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Abstract. I discuss the importance of dynamical models for exploiting survey data, focusing on the advantages 
of "torus" models. I summarize a number of applications of these models to the study of the Milky Way, including 
the determination of the peculiar Solar velocity and investigation of the Hyades moving group. 



1 Introduction 

Current studies of the structure of the Milky Way are dom- 
inated by a series of major observational programs running 
from ESA's Hipparcos mission |1 1, through major photo- 
metric surveys such as the SDSS |2 1, spectroscopic surveys 
such as RAVE |3|, to ESA's scheduled Gaia mission ||4], 
which aims to return photometric and astrometric data for 
10'' stars and low-dispersion spectra for ~ 10^ stars. 

Turning these data-sets into a consistent picture of the 
current structure and the assembly of the whole Galaxy, in- 
cluding the dark-matter content, is an ambitious and impor- 
tant goal. It is likely to be impossible without sufficiently 
sophisticated models that can be used to interpret the data 
(for example compensating for the observational biases of 
the various surveys). 

Models of the gross structure of the Galaxy have been 
produced with varying levels of complexity. Mass mod- 
els |5 6| simply give the density distribution of the various 
components of the Galaxy, and thus the Galactic potential. 
Kinematic models, such as those produced by galaxia Q, 
specify the density and velocity distributions of the lumi- 
nous components of the Galaxy, but do not consider the 
question of whether these are consistent with a steady state 
in any Galactic potential. The Besangon Galaxy model lO 
is primarily a kinematic model (in that it is not constrained 
by Newton's laws of motion on large scales), with a dy- 
namical element used to determine the vertical structure of 
the disc. Fully dynamical models |9| have a distribution 
of stars in phase-space which is made from phase-mixed 
orbits in a given Galactic potential, and therefore a distri- 
bution function (df) which is a function of the integrals 
of motion, which by Jeans theorem ifTOl means it is in a 
steady-state. 



2 Benefits of dynamical models 

A previous paper [9| provides a detailed discussion of the 
benefits of dynamical models. Here I just sum up the main 
points. 
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The most obvious advantage of dynamical models is 
that, unlike kinematic models, they allow us to deduce the 
gravitational potential of the Galaxy. Existing mass mod- 
els were fit to observations under the assumption of near 
circular orbits for various tracers, which is only suitable 
for a small subset of astrophysical objects, found close to 
the Galactic plane. Exploiting richer data-sets requires far 
more careful modelling. This will allow us to infer the dis- 
tribution of dark matter in the Galaxy, under the assump- 
tion that the Galaxy is approximately in a steady state. 

A second, complementary, advantage is that dynamical 
models allow us to connect objects that we can observe to 
objects those we can not. This means we can use observa- 
tions in the solar neighbourhood to learn about the struc- 
ture of the Galaxy at large. For example IfTTII showed that 
if the stellar halo were in virial equilibrium, more than half 
the stars of the stellar halo would be on orbits that bring 
them through the solar neighbourhood. 

An additional advantage of dynamical models is that 
the associated df depends only on the three integrals of 
motion (either explicitly or implicitly), as opposed to the 
full six dimensions of phase-space. 



3 Torus models 

Dynamical modelling has been dominated by Schwarzschild 
modelling [I2|, especially for analysing the dynamics of 
early-type galaxies. This technique involves first integrat- 
ing a number of orbits in the adopted gravitational poten- 
tial and then seeking weights for these orbits such that 
the weighted sum of the orbits reproduces the observa- 
tional data. More recently, the "made-to-measure" (M2M) 
technique introduced by LISJ has been used to produce A^- 
body models which can be fitted in a broadly similar way 
II14I1 51. though in this case the particle weights (which are 
eff'ectively the orbit weights) are determined "on-the-fly", 
rather than after the orbit has been integrated. 

Both Schwarzschild and M2M models have dfs which 
are implicitly function of the integrals of motion, as they 
are constructed from phase-averaged orbits. An alternative 
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Fig. 1. Distribution of stars in vn (left) and (right) in tiie Solar neighbourhood, data from the Geneva-Copenhagen survey assuming 
the value of Vq given in eq.|2](histogram), and model from 1 17] which is a single df /(J) which must fit both the vj^ and distributions 
(curves, solid curve is from the full model, dashed curve is thin disk contribution only). The dotted vertical line in the plot is the 
assumed circular speed at the Sun. Clearly the distributions are significantly offset from one another. The only way to bring the data 
and model to reasonable agreement is to apply a correction of ~ 7kms"' to Vq. Figure reproduced from |T7| with permission. 



modelling strategy is to produce models which are explic- 
itly functions of the integrals of motion. It is not possible 
to produce plausible models of the Galactic disc with dps 
which are functions of the classical integrals, f{E, L^), as 
these models have equal velocity dispersions in the radial 
and vertical directions, which is not the case in the So- 
lar neighbourhood 1 16|. However, one can instead use the 
three orbital "actions", 7,, as the integrals of motion and 
use analytic dps, /(J), which are appropriate for the Galac- 
tic disc (e.g. IflTl ). 

The three actions 7, and three conjugate angles coordi- 
nates Oi provide canonical coordinates for six-dimensional 
phase space IfTSll . The conventional phase space coordi- 
nates w = (x, v) are 27r-periodic in the angles. The actions 
are conserved quantities for any orbit, and the angles in- 
crease linearly with time: 

e(t) = 6(0) + Q(J)t, (1) 

where the components of £2 are the orbital frequencies. 

The major obstacle to using dps of the form /(J) is 
that the relationship between phase space coordinates w 
and J, 9 is only known analytically for a very limited set 
of gravitational potentials, none of which provides a rea- 
sonable approximation to the Galactic potential. There are 
two available approaches for determining the relationship 
between w and J: 

- The adiabatic approximation, in which motion in the 
radial and vertical directions are largely decoupled 1171191 . 

- Torus modelling, in which the relationship between J, 6 
and w in an isochrone potential (for which it is known 
analytically [18J) is numerically distorted to fit the po- 
tential of interest 0201211221 . 

Comparison of these two approaches has shown that for 
most purposes they agree to reasonable accuracy up to as 
far as ~ 2.5 kpc from the Galactic plane 1 19|. The adiabatic 
approximation has the advantage that it does not require 



specialised torus-fitting computer code, and can straight- 
forwardly determine the value of J for a given w. Torus 
modelling finds all of the values of w associated with a 
given J, but can only find J given w as an iterative process. 
Torus modelling has the advantages that it can tell us about 
the coupling between different components of motion (e.g. 
the tilt of the velocity ellipsoid [19.1 ), and that it allows us 
to find the angle variables 0. 

4 The local standard of rest 

As mentioned in Section |2l one major advantage of using 
dynamical models is that the dimensionality of the models 
is reduced by the assumption that the Galaxy is made up of 
phase mixed orbits - /(J) depends on only three actions as 
opposed to the six dimensions of /(w). The value of this 
simplification of the model is well illustrated by a recent 
revision in the peculiar Solar velocity relative to the local 
standard of rest, Vq. 

The value of v© was found by |23| using observations 
by Hipparcos. The components of velocity towards the Galac- 
tic centre (U), in the direction of Galactic rotation (V) and 
towards the north Galactic pole (W) were analysed sepa- 
rately. The V-component of the Solar velocity is the most 
difficult to determine as asymmetric drift IfTSlI means that 
average stellar velocity lags the circular velocity. Stromberg's 
equation was used by |23 1 to extrapolate from the observed 
populations (separated by colour) to a hypothetical popula- 
tion with zero velocity dispersion, which would have zero 
asymmetric drift. The value of Vq found, 

Uq, Vq, Wq = (10.00+0.36, 5.25+0.62, 7.17±0.38)km s^^ 

(2) 

was the widely accepted value for over a decade. 

Figure [T] shows the distribution of stars in and ur 
in the Solar neighbourhood, determined from Hipparcos 
observations and the Geneva-Copenhagen survey Il24ll as- 
suming this value of Vq (histogram), and the best fitting 
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Fig. 2. Frequencies & (left) and a convenient projec- 
tion in action space (right) for stars in the solar neighbourhood 
from a simulation of the accretion of a satellite galaxy, taken 
t = 1.5 Gyr after the satellite was accreted. The stars form patches 
in frequency space with spacing AQ^ between patches such that 
AQ^t/2n ~ 1 (the spacing in i3,- is not as simple because a 
two ranges of 9^ correspond to the solar neighbourhood). The 
same separation into discrete clumps is visible in action because 
i2 = These figures are taken from |22|. 

analytic df /(J) from ifTTl (solid curve). This gives new 
insight because the distributions in U and V are not con- 
sidered as if they were independent, but instead it is recog- 
nised that a single dynamically consistent df must fit both. 
The only physically plausible way to bring the model and 
data distributions into agreement is to apply a correc- 
tion of ~ 7kms ' to the value of (altering the circular 
velocity moves both distributions and has a negligible ef- 
fect). This is ~ 11 cr from the widely accepted value! This 
contention, that the previously assumed value of Vq must 
be in error, has since been supported by analysis of Galac- 
tic maser sources with measured parallaxes |25|, and ex- 
plained as an error in the application of Stromberg's equa- 
tion, associated with the metallicity gradient in the disc 



5 Beyond steady-state models 

Axisymmetric models with dfs of the form /(J) cannot 
fully describe the Galaxy, as it is not in fact in a steady 
state, but they are an important first step towards interpret- 
ing observations of our Galaxy, and it is likely to be very 
fruitful to study observational data for structures that can- 
not be explained by these models. We can then look for 
explanations of these structures as signatures of other fea- 
tures, such as the Galactic bar, spiral or warp, or matter that 
has been accreted. Some examples of using torus models 
to explore these signatures already exist. 



5.1 Signatures of accreted satellites 

The appearance in angle-action coordinates of an accreted 
satellite galaxy was explored by |22|. Long after phase 
mixing has rendered an accreted satellite indistinguishable 
from the background population in position, there is a strong 
relationship between the stars' positions and their orbital 
frequencies (because all the stars were once collected in 



the same small volume, when they were part of the satel- 
lite). This means that a sample of these stars taken from 
a finite volume is only found in certain small volumes in 
frequency space. 

In Fig|2]I show figures from f22] of the frequencies and 
actions of stars in a finite volume about the solar position 
in a simulation of the disruption of a satellite galaxy. In 
both cases the figure shows the distribution 7.5 Gyr after 
the satellite was disrupted. The distribution in frequency is 
clearly divided into finite "patches", and the distribution in 
action is even more cleanly divided, because £2 = £2{J) and 
it provides a more convenient projection of the distribution. 

By considering the spacing between the "patches" in 
frequency space, along with the angles of the individual 
stars, 1 22 1 showed that it was even possible to determine 
with high accuracy the time since the satellite was dis- 
rupted. Similar techniques could even be used to determine 
the potential of the Galaxy (as the potential must allow the 
stars to all have come from the same initial satellite). Simi- 
lar work has been caiTied out by other authors showing that 
using the orbital frequencies alone one can achieve some 
of these objectives, even for cosmological simulations with 
numerous accreting satellites and a non-static background 
potential IIZTll . 



5.2 Signatures of Lindblad resonances 

A recent study |28 1 showed that in addition to having a re- 
lationship of the form IQr + mQ^ ^ const, stars that have 
recently been trapped at a resonance with a perturber also 
follow a relationship of the form 10,- + mO^ ^ const, where 
in both cases I and m are integers, with the perturbation 
having m-fold symmetry and / being - 1 for an inner Lind- 
blad resonance and +1 for an outer Lindblad resonance. 
The Hyades moving group, which is a very strong feature 
of the local velocity distribution |29|, lies around a straight 
line in the J^, Jr plane (Fig [3]), of the sort associated with 
the condition IQ,- + mQ^ ^ const (a resonance line in action 
space). This finding is indeed consistent with a resonance 
line for the both the cases / = +1 and a range of values 
of m (with the exact details of the resonance lines depend- 
ing quite sensitively on the details of the Galactic potential 
assumed). 

It was claimed by ||28l that the distribution of stars in 
angle coordinates clearly indicated that the Hyades were 
associated with an inner Lindblad resonance. This is be- 
cause, in angle space, the stars of the Hyades moving group 
are associated with an overdensity in the quantity -Oy+mO^ 
for various values of m, and this overdensity did not appear 
to significantly shift in position as a function of Jr. 

In 1301 1 compared the observed structure to torus mod- 
els, and demonstrated the significant and non-intuitive im- 
pact of selection effects. In Fig[3]I show the distribution of 
Solar neighbourhood stars in the 0,-, 0^ plane. Selection ef- 
fects are responsible for the overall structure of the density 
distribution, most notably the high densities around 0^ - 0, 
0r - or +7r and the near absence of stars with 0r < 0, 
0^ > or 0r > 0, 00 < 0. A less obvious selection eff'ect 
brought to Ught by torus models is that stars with a given 
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(with resonant conditions on both J and 0) indicated that 
the Hyades overdensity was consistent with stars trapped 
at either an inner or outer Lindblad resonance. 
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Fig. 3. Distribution of stars in action (upper) and angle (lower) in 
the Solar neighbourhood with positions and velocities given by 
the Geneva Copenhagen survey. The gross structures of the dis- 
tributions are due to selection effects. The Hyades moving group 
can be seen as an overdensity in action that is spread out in 7, at 
around = 0.97, tending towards slightly lower with increas- 
ing Jr (consistent with a resonance line in ii), and as an over- 
density in angle around 9r = -n/2. The expected angle relation 
for stars at resonance would produce overdensities that tend to 
increased 6, with increased 6^ for inner Lindblad resonances, and 
to decreased 9^ with increased 6^ for outer Lindblad resonances. 
However, selection effects reshape overdensities in angle space 
quite significantly. Figures are adapted from 1301 



value of J, found in the Solar neighbourhood are also very 
strongly restricted in their possible range of 0. 

Using these models I was able to show that observed 
overdensities in 6r+m0^ (which should coiTespond to outer 
Lindblad resonances) were not the result of selection ef- 
fects (as claimed by ||28l ). Also, these selection effects 
mean that the stars associated with the resonance line in 
action space have to lie near certain lines in angle space, 
otherwise they will not be observed in the Solar neighbour- 
hood. It is the interplay of the resonant conditions on J and 
6, and the fact that the stars lie in the Solar neighbour- 
hood that determines the distribution in both angle and 
action, and considering either distribution independently 
of the other can lead to false conclusions. The examina- 
tion of torus models which included resonant components 
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